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摘 要 : 针对 基于 Gabor 小 波幅 值 与 相位 的 人 脸 特 征 提取 方法 的 特征 级 联 方式 使 得 特征 向 量 维度 较 高 的 问题 ， 提 出 
了 一 种 改进 的 Gabor 小 波 变换 特征 提取 算法 。 该 算法 计算 局 部 幅 值 特征 和 局 部 相位 特征 ， 增强 了 每 个 像素 的 局 部 关 

联 性 ; 然后 通过 实验 选 定 加 权 系 数 ， 将 幅 值 特征 与 相位 特征 进行 加 权 融 合 。 实 验 结果 表明 ， 所 提出 的 算法 与 改进 前 

的 算法 相 比 ， 降 低 了 特征 向 量 的 维度 ， 且 提高 了 最 终 的 人 脸 识别 率 。 
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Feature extraction algorithm based on improved Gabor wavelet transform 
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College of Computer, Nanjing University of Posts & Telecommunications, Nanjing 210003, China) 


Abstract: Aiming at the problem that the feature cascading method for face feature extraction based on Gabor wavelet 
amplitude and phase makes the feature vector dimension higher, this paper proposed an improved Gabor wavelet transform 
feature extraction algorithm. The algorithm calculated local amplitude and local phase features, and enhanced the local 
correlation of each pixel. Then the algorithm weighted fusion for amplitude and phase features by weighting coefficients 
which experiment selected. Experimental results show that the proposed algorithm reduces the dimension of the feature 
vector and improves the final face recognition rate compared with the pre improved algorithm. 
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Gabor 小 波 变换 特征 提取 算法 。 该 算法 计算 每 个 像素 点 的 局 


0 Sla 部 幅 值 特征 和 局 部 相位 特征 来 分 别 代替 该 像素 点 的 幅 值 特征 
特征 提取 是 人 脸 识 别 系统 中 最 关键 的 技术 之 一 ， 指 的 是 。 和 相位 特征 ， 增 强 每 个 像素 点 的 局 部 关联 性 ， 采 用 加 权 融合 

通过 计算 分 析 人 脸 图 像 得 到 具有 一 定 区 分 度 的 数据 来 描述 人 的 方式 对 幅 值 特征 与 相位 特征 进行 融合 ， 降 低 人 脸 特征 向 量 

从 特征 的 过 程 。 但 是 目前 已 有 的 人 脸 特征 提取 方法 在 算法 对 ”的 维度 并 提高 人 脸 特征 的 表征 能 

于 光照 与 姿态 等 外 界 因素 干扰 的 情况 下 ， 算 法 提取 的 特征 的 

名 标 性 不 高 ， 还 不 能 适用 于 一 般 情况 。 所 以 ， 如 何 找到 一 种 。 1 ”相关 工作 

RISEREU EHE EA EDERE JEA SCE ERR 目前 基于 Gabor 小 波 变换 的 人 脸 特征 提取 算法 从 图 像 变 

问题 换 区 域 大 小 上 可 以 分 为 解析 法 和 整体 法 中。 
0 解析 法 即 计算 每 一 个 特定 标记 区 域 的 Gabor 小 波 特征 ， 

旬 越 来 越 多 ， 受 到 包括 光照 、 姿 态 、 庶 挡 等 外 界 因素 的 影响 。 所 以 该 类 方法 需要 选 定 一 些 特征 域 或 者 特征 点 (如 眼睛 、 自 


壁 如 当 光 照发 生 较 大 变化 时 ， 同 一 个 人 的 人 脸 图 像 可 能 比 在 T 嘴巴 等 ), 然后 求 得 各 个 标记 区 域 的 Gabor 特征 并 将 其 组 
相同 光照 条 件 下 的 不 同人 的 人 脸 图 像 的 变化 大 很 多 ， 因 为 光 合成 “Gabor jet". Lades 等 人 中 最 早 提出 Gabor 小 波动 态 连 
照 的 变化 使 得 处 于 光照 较 暗 下 的 人 脸 的 纹理 信息 被 遮盖 ， 所 ， 接 结 构 并 将 其 应 用 于 人 脸 识 别 中 ， 该 方法 将 人 脸 表示 为 格 状 
以 导致 同一 个 人 的 人 脸 特征 发 生 较 大 差异 。 当 人 脸 姿 态 变化 HARAR, MERMER Gabor 小 波 变换 得 到 的 特征 向 量 标 
局 度 过 大 或 者 人 脸 部 分 信息 被 谈 挡 时 ， 会 导致 人 脸 部 分 关键 。” 记 图 的 节点 ， 用 连接 节点 的 距离 向 量 标 记 图 的 边 ， 根 据 图 像 
特征 缺失 ， 从 而 无 法 提取 完整 的 人 脸 特 征 。 因 此 ， 传 统 的 特 。 中 的 节点 与 连接 之 间 的 相似 性 进行 人 脸 识 别 ， 但 是 在 图 像 中 
征 提 取 算 法 对 于 此 类 图 像 的 特征 提取 的 效果 不 够 理想 。 而 通 ”搜索 特征 点 的 计算 量 比 较 大 。 为 了 解决 这 一 问题 ，Wiskott 
过 研究 发 现 Gabor 小 波 类 似 于 哺乳 动物 视觉 系统 中 视觉 神 ”等 人 中 使 用 弹性 图 匹配 方法 对 人 脸 的 关键 部 位 进行 标定 ， 然 
经 元 ， 在 局 部 空间 内 有 很 好 地 提取 信息 的 性 能 ， 由 于 Gabor ”后 针对 关键 部 位 进行 了 Gabor 小 波 变换 ， 最 后 将 多 个 关键 部 
小 波 具 有 良好 的 方向 选择 性 、 多 尺度 性 以 及 对 于 光照 、 姿 态 。 位 的 Gabor 小 波 特征 级 联 。 在 此 过 程 中 弹性 图 的 匹配 较为 耗 
等 变化 的 不 敏感 性 ， 使 得 Gabor 小 波 在 图 像 处 理 和 模式 识别 ”时 。 因 此， 为 了 避免 弹性 图 的 构建 ，Kalocsai SAS AS 
中 得 到 了 广泛 的 应 用 请 和 。 的 关键 特征 点 直接 进行 Gabor 小 波 变 换 ， 再 将 这 些 特征 级 联 
因此 ， 为 了 进一步 提升 基于 Gabor 小 波 变换 特征 提取 算 ”起 来 .但 是 解析 法 都 需要 对 人 脸 的 关键 特征 点 进行 手工 标定 ， 
法 提取 特征 的 表征 能 力 ， 并 降低 特征 维度 ， 本 文 提出 改进 的 其 过 程 较为 繁琐 。 
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整体 法 则 是 将 图 像 整体 进行 Gabor 小 波 变换 ， 不 再 需要 F (@)=1(2)*G,, (2) 3) 
对 特征 点 进行 手工 标记 。Liu 等 人 外 使 用 多 尺度 与 多 方向 的 M 
Gabor 核 函 数 对 完整 的 人 脸 图 像 进行 卷 积 ， 将 提取 的 Gabor 加 E E 
其 中 : E.G) 表示 bor 小 ; ESI 和 特征 图 像 ， 1CO 表示 
幅 值 特征 作为 人 脸 的 特征 向 量 。Shan 等 人 no 将 Gabor 幅 值 特 T GREOPQSPUSHPUB SINT ESDEGS EIS USUS 
征 与 Boosting 策略 相 结合 ， 降 低 特征 维度 得 到 AdaGabor 特 ”初始 图 像 ，* 表示 卷 积 运算 。 
征 ， 最 终 通过 此 特征 对 人 脸 进行 识别 。 此 前 的 这 些 方法 提取 
的 都 是 Gabor 的 幅 值 特征 ， 忽 略 了 Gabor 的 相位 特征 。 事 实 EM. Lu. N 
E, Gabor 相位 特征 中 含有 丰富 的 可 以 用 于 分 类 识别 的 信息 ， 
经 过 研究 表明 ，Gabor 相位 特征 对 于 光照 等 外 界 因 素 的 变 
比 具有 较 高 的 鲁 棒 性 中 1。Zhang 等 人 0 提出 了 基于 Gabor 小 10| 5|20|20|20 -10| 5 |15| 0 | 0 
波 相位 信息 的 直方 图 模型 ， 该 方法 使 用 了 Gabor 小 波 变换 、 NM E NEC UE 
空间 区 域 直 方 图 等 ， 使 得 该 方法 对 于 外 界 因素 的 变化 具有 良 ————Ma emm 
好 的 鲁 棒 性 。 为 了 更 好 地 提升 基于 Gabor 小 波 变换 的 特征 提 GESPSERES aa ERE 
取 算 法 的 鲁 棒 性 ， 文 献 [12] 提 取 了 人 脸 的 Gabor 幅 值 特征 与 
相位 特征 ， 并 将 两 者 级 联 获得 最 终 人 脸 特 征 向 量 。 但 是 目前 ea ee a ES o 
基于 Gabor 小 波幅 值 和 相位 信息 的 特征 提取 方法 一 方面 较 少 图 1 二 维 卷 积 运算 
考虑 幅 值 特征 的 局 部 邻 域 之 间 的 关系 ， 使 得 幅 值 特征 的 局 部 Fig.1 Two dimensional convolution operation 
关联 性 较 弱 ， 另 一 方面 合成 幅 值 特征 与 相位 特征 的 方式 大 多 、 NN 
为 串联 ， 所 以 使 得 最 终 的 人 胎 特 征 向 量 的 维度 较 高 且 表 征 能 了” 改进 的 Gabor 小 波 变换 特征 提取 算法 
力 不 能 达到 最 好 。 3.1 Gabor 局 部 幅 值 特征 
综 上 所 述 ， 解 析 法 由 于 需要 手工 标定 特征 ， 所 以 其 过 程 基于 Gabor 幅 值 统计 特性 的 特征 提取 算法 所 提取 的 图 像 
较为 繁杂 。 整 体 法 中 基于 Gabor 小 波幅 值 信息 的 特征 提取 方 。” 幅 值 特征 信息 相 较 于 图 像 全 局 特征 受 外 界 因素 的 干扰 要 弱 一 
法 对 于 限定 条 件 下 的 人 脸 特征 提取 有 较 好 的 效果 ， 因 为 其 通 。 " Meu 
T 4E Gabor h C CO 为 复数 形式 ， 所 以 其 实 当 
过 滤波 得 到 的 幅 值 信 息 可 以 很 好 地 反映 人 脸 的 纹理 特征 ， 然 a Sor cde 和 
而 复杂 条 件 下 的 人 脸 局 部 的 纹理 特征 可 能 会 发 生 改变 ， 所 以 p 
e(G,,(z)) 和 虚 部 Im(G, JG 分 5 p 
导致 基于 Gabor 小 波幅 值 信息 的 特征 提取 方法 对 于 受 外 界 环 和 NURSES 
境 干 扰 的 人 脸 提取 的 特征 的 表征 能 力 不 强 。 而 基于 Gabor 小 - 
波幅 值 和 相位 信息 的 特征 提取 方法 融入 了 对 于 光照 等 外 界 因 Re(0 o) Ste ( e a-ev|-T] ^ 
素 具有 较 强 鲁 棒 性 的 相位 信息 ， 增 强 了 特征 的 表征 能 力 ， 但 
是 两 种 信息 的 级 联 方式 导致 特征 维度 较 高 。 
" G, (Zz Sarex sin(k, 
2  —H Gabor 小 波 变换 In(6,,(9))- (x a Jen C 
Gabor 小 波 是 由 Gabor 函数 经 过 尺度 伸缩 和 旋转 生成 的 " i 
i 第 2 章 中 式 G), 可 以 得 到 图 像 的 Gabor 特征 Fo (2) , 
一 族 复 函数 系 ， 具 有 很 好 的 视频 局 部 化 特征 和 多 分 辩 性 能 ， 人 
能 够 提取 图 像 局 部 细微 变化 031。 对 于 图 像 进 行 Gabor 滤波 ， me ————á 
SERRE Gabor 小 波 变 换 。 二 维 Gabor aEBunmasug IU Ie CO 也 为 复数 ， 其 实 部 为 Re(Re GO ENDS E, D 
DHrZRUM, z 
n 比 幅 值 特征 为 M), ROLLS 
k, kis : o 
a. ()- ren [- 2o Jes. (- 2) (D M, (2) = (Re CE, (Y UmCE,, (2) (6) 
Kd oz 代表 图 像 中 某 个 像素 点 在 空间 中 的 坐标 值 ， 可 以 用 通过 卷 积 运算 获得 图 像 的 Gabor 幅 值 特征 后 ， 为 了 增强 
B "NEM 局 值 特征 的 局 部 关联 性 ， 则 计算 每 个 像素 点 的 局 部 幅 值 特征 
= 不 ; KE... 示 滤 波 器 的 中 心 频率 。 S en à 
EYRE: Kelkan] 表示 滤波 器 的 中 心 频率 来 代替 该 像素 点 的 幅 值 特征 。 本 文选 取 3x 3 大 小 的 窗口 ， 作 
" H e 为 一 个 像素 点 的 邻 域 ， 计 算 该 窗口 中 周围 八 个 像素 点 的 平均 
L9 275 -uc, HH 20,1,2,..., 不 bor J&i E 
korti wor. 0,12…7) 表 不 Gabor WA — 幅 值 ， 用 该 平均 幅 值 来 代 蕉 窗口 中 心 像素 点 的 幅 值 。 当 像素 
的 方向 vv=0,1.2,3, 和 表示 尺度 ; k 表示 方向 的 总 数目 。i 点 位 于 幅 值 特征 图 谱 的 边界 时 ， 邻 域 不 在 幅 值 特征 图 谱 范 轩 
是 复数 运算 符 ，* 表示 求 对 应 值 的 欧 氏 范 数 ，s 决 定 了 高 斯 ” 内 的 幅 值 补 零 。 假 设 像素 点  ， 其 邻 域内 各 像素 点 
窗 宽 与 波长 之 比 ， 通 常 e=2x 时 实验 效果 最 佳 。 
'NDEUL Gubor 特征 图 谱 可 以 通过 图 像 避 Gabor ak TETO O AMR 2 所 示 ， 则 其 幅 值 Mi,() 的 计算 公式 
的 核 函 数 进行 卷 积 运算 而 获得 。 如 图 1 所 示 ，A 可 以 看 做 待 、 m T 
7 VAR 有 Hur fü 8 H 特征 图 X M, (z o 
处 理 的 图 像 ，B 可 以 看 做 Gabor 小 波 的 核 函数 。 卷 积 运算 的 如 式 (7) 所 示 ， 最 后 计算 获得 局 部 幅 值 特征 图 谱 M, (z) 
ANS M (z 2M. (2) (7) 
C(x.») -3YA(Q. ) Bx-i.y - j) (2) 8 
i 为 了 能 够 有 效 地 表征 图 像 全 局 特征 以 及 降低 特征 维度 ， 
那么 ， 将 输入 的 人 脸 图 像 进行 Gabor 小 波 变换 ， 则 运算 。 ”将 得 到 的 局 部 幅 值 特征 图 谱 进 一 步 划 分 为 若干 个 矩形 不 重 车 
可 以 定义 为 且 大 小 相等 的 子 块 ， 子 块 数 记 为 K ， 然 后 分 别 计算 各 子 块 的 
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统计 特征 。 统 计 特征 选择 样本 均值 和 标准 差 这 两 个 最 为 典型 
的 统计 量 。 
Z1 Z2 Z3 
Zg Zo Z4 
Z7 Z6 Zs 
图 2 邻 域内 各 像素 点 分 布 示意 图 


Fig.2 Schematic diagram of distribution of each pixel in 


neighborhood 


假设 一 幅 图 像 工 大 小 为 到 xm ， 
图 谱 划 分 了 个 大 小 相同 的 


针对 Gabor 局 部 幅 值 特征 
子 块 ， 每 个 子 块 图 像 为 


1, (kK-0,1,..., K -1)， 每 个 子 块 图 像 内 的 均值 M4 和 标准 差 


Mia 的 计算 公式 分 别 为 


(8) 


(9) 


则 图 像 I 的 幅 值 特 征 可 以 表示 为 


H. rH 


3.2 Gabor 局 部 相位 特征 
于 Gabor 小 波 的 相位 特 各 


((MZ, Mz, )|u 9 0.1... U, v 20.1... VK 20,1... K —1} 


(10) 


tr: U 和 V 分 别 表 示 Gabor 小 波 核 函 数 的 方向 和 尺度 数 。 


E 对 于 光照 有 着 很 强 的 分 辨 


性 ， 所 以 本 章 根据 Daugman 方法 05 提 取 图 像 经 过 Gabor 小 


波 变换 后 的 局 部 相位 特征 ， 将 其 与 局 部 幅 值 特征 进行 融合 ， 


得 到 对 于 光照 具 
Daugman 方法 将 经 过 二 维 


可 更 强 鲁 棒 性 的 特征 。 
Gabor 小 波 变换 后 的 图 像 中 的 


每 一 个 像素 点 的 信息 进行 二 进 制 转 换 ， 即 根据 每 个 像素 点 的 


实 部 信息 和 虚 部 信息 的 正 负 来 确定 其 二 进 制 的 值 ， 转 换 规 则 


如 式 (1D 和 式 (12) 所 示 05。 


zo 


as | 


0, Re(F,, (z))»0 
l, Re(F,,(z)) «0 


0. Im(F,,(z))>0 
L Im(F,,(z)) «0 


(41) 


(12) 


其 中 : z 表示 图 像 中 某 个 像素 点 化 标 ;Re(,(z)) 、Im(,(z)) 


分 别 表 示 像 素 点 的 实 部 信息 和 虚 部 信息 ， P). Pie (2) 分别 


表示 像素 点 的 实 部 信息 和 虚 部 信息 的 二 进 制 转换 值 。 图 像 中 
每 一 个 像素 点 的 实 部 信息 和 虚 部 信息 组 合成 的 坐标 


(Re(R (z)).Im(F., (z))) 都 可 以 在 复 
图 3 所 示 。 


图 中 ， 


6, G) ET 28 T BRIV 象限 当 Re( 


坐标 系 内 找到 对 应 的 点 , 如 


90,,(z) 表示 华 标的 相位 角 。 当 Re(F,, (2)) 2 0 HF, 


F) SORT, Po (7) 位 于 第 
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IL EI 象限 ; 同 理 可 以 推出 Im(F, (2)) M 6,. (2) 的 关系 。 所 以 
式 (11) 和 (12) 可 以 分 别 变 形 为 式 (13) 和 (14)。 


«La [O Qs)e {LIV) 
pA te em 
«cS [O 8. etm 
An "y 


式 (13) 和 (14) 只 是 分 别 对 实 部 与 虚 部 信息 进行 了 二 进 制 
转换 , 并 没有 对 实 部 与 虚 部 信息 之 间 的 关系 进行 二 进 制 转 换 。 
为 为 在 复 平 面 中 , 相位 角 6. G) RI Re(F, (2) « Im(F, (2) 存在 
式 (15) 所 描述 的 关系 。 

m T 

FA, BXO3 OA) RE, EI). Ph (z) 相等 ， 即 


9,(7) 位 于 第 I 或 王 象 限时 , tan(6,,(2)) 0; 24 E (z) Rel) 


不 相等 ， 即 9,,(z) 位 于 第 I 开 或 IV 象限 时 ，tan(6,,(z))<0。 所 


以 按照 同样 的 三 进 制 转换 规则 对 tan (8, (2)) 进行 二 进 制 编码 


得 到 Pw(z) ， 如 式 (16) 所 示 。 


Jo .Ge 下 
C NPIPIA MS 
T Pr (2) = PË (2) XORB/" (2) a7) 


为 了 增强 每 个 像素 点 与 其 局 部 邻 域内 像素 点 的 相位 特征 
的 关联 性 ,同样 选 定 3x3 大 小 的 窗口 ， 作 为 一 个 像素 点 的 邻 
域 。 将 窗口 内 其 他 八 个 位 置 的 二 进 制 编码 与 中 心 像素 点 的 编 
码 进行 异 或 运算 得 到 新 的 八 个 二 进 制 编码 ， 然 后 将 这 八 个 二 
进 制 编码 按照 顺序 排列 成 一 个 8 位 二 进 制 数 ， 最 后 将 其 转换 
为 一 个 十 进 制 数 来 表示 中 心 像 素 点 的 相位 特征 。 当 像素 点 位 
于 相位 特征 图 谱 的 边界 时 ， 邻 域 不 在 相位 特征 图 谱 范围 内 的 
相位 特征 编码 补 零 。 具 体 运算 过 程 如 图 4 所 示 。 


A 
I —210 1—200 


(Re(F, (2) Im(F, (2) 


M—>l1 IV—501 


图 3 REMATE R 


Fig.3 Schematic diagram of complex coordinate system 


UN Zo 
TEDA oloi 0f112 
ol YI1i XR. mm Po 1| |3 
ucc "^^^ ool ^ ^^— aba 
二 进 制 编码 : 00101101 了 Puv(Zo)=45 


图 4 局 部 相位 特征 运算 过 程 


Fig.4 Operation process of local phase feature 
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RP, 为 某 一 像素 点 坐标 ， 其 邻 域 内 各 像素 点 


zi(i=1,2,.…,8) ， 像 素 点 z 的 局 部 相位 特征 为 , (zo) ， 最 终 得 


到 图 像 的 局 部 相位 特征 图 谱 Aa (2) 。 


与 3.1 节 的 局 部 幅 值 特征 提取 类 似 ， 为 了 有 效 地 表征 图 
像 全 局 特征 以 及 降低 特征 维度 ， 一 幅 大 小 为 mxn 的 图 像 I 通 


过 二 维 Gabor 小 波 变换 得 到 局 部 相位 特征 图 谱 R.(z) 后 ， 将 


其 进一步 划分 为 若干 个 矩形 不 重叠 且 大 小 相等 的 子 块 ， 子 块 
数 记 为 Kk， 然后 分 别 通过 式 (18) 和 (19) 计 算 各 子 块 


1, (k-0,1,..., K -1) 的 均值 Ph 和 标准 差 Prue o 


P, (z 
"uo (18) 

(mxn)/K 
pe, |e POR) ig; 

i (mxn)/ K 

则 图 像 T 的 相位 特征 可 以 表示 为 
(25, 05,)|u201...U.v201.., Vk-01...K-H] (20) 


其 中 : UM V 分 别 表示 Gabor 小 波 核 函数 的 方向 和 尺度 数 。 
3.8 ”特征 融合 

特征 融合 思想 是 对 多 种 特征 进行 综合 分 析 、 评 判 和 决策 
处 理 的 过 程 ， 即 通过 融合 得 到 一 个 使 各 特征 对 结果 影响 贡献 
最 大 的 综合 特征 。 目 前 对 于 Gabor 小 波 变换 的 幅 值 特征 与 相 
位 特征 的 融合 最 常 采用 的 方式 为 级 联 ， 即 串联 法 ， 这 种 融合 
方法 是 将 各 个 特征 视 为 一 种 假设 的 平等 ， 认 为 每 一 个 特征 对 
于 最 终 的 识别 结果 具有 同等 的 影响 ， 这 种 方法 并 不 能 使 各 个 
特征 对 识别 正确 的 贡献 作用 发 挥 到 最 大 。 而 特征 加 权 融 合 方 
法 是 为 每 一 个 特征 定义 权 值 ， 通 过 实验 来 确定 各 个 特征 的 贡 
献 值 。 令 幅 值 特征 的 权 值 为 W ， 相 位 特征 的 权 值 为 W, ， 设 
定 每 一 个 特征 的 权 值 范围 为 01~0.9， 且 W,+W=1。 所 以 最 
终 的 融合 特征 计算 公式 如 式 (21) 所 示 。 


MIX „vx =W, (Mi ? My. ) +W, (Pe; , Pam ) (2 1) 


u 


二 


中 : u=0,1,..,U; v=0,1,..,V ; k=0,1,...,K-1; U, VĀKS 
别 表 示 Gabor 小 波 核 函 数 的 方向 和 尺度 数 以 及 特征 图 谱 的 分 
块 大 小 。 
3.4 算法 描述 
改进 的 Gabor 小 波 变换 特征 提取 算法 主要 由 四 部 分 组 
成 ， 首 先 需要 对 人 脸 图 像 进行 二 维 Gabor 小 波 变换 ， 其 次 提 
取 每 个 像素 点 的 幅 值 特征 并 计算 每 个 像素 点 局 部 幅 值 特征 ， 
然后 提取 每 个 像素 点 的 相位 特征 并 计算 每 个 像素 点 局 部 相位 
特征 ， 最 后 对 局 部 幅 值 特征 和 局 部 相位 特征 进行 加 权 融 合 。 
改进 的 Gabor 小 波 变换 特征 提取 算法 的 伪 代 码 如 算法 1 所 示 。 
算法 1 改进 的 Gabor 小 波 变换 特征 提取 算法 
输入 一 幅 大 小 为 mxn 的 人 脸 图 像 ， 输 出 人 脸 图 像 的 特征 
向 量 ; 
l. za 擂 第 一 个 像素 点 的 坐标 ; 


2.fori 人 lto mxn; 


3.  F.(z)€1(z)*G, (a); 


4. Re(F,, (z)) € I (z)*Re(G,, (z)) $ 
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5. — Im(E, (z)) € 1(z)*Im(G,, (5); 


储 小 波 变换 


M 


6. end for// 对 人 脸 图 像 进行 
7.fori*-1to mxn; 


8. — M, (x) - A(Re(F, (5) Qm(E, (5). : 
9. forjc€lto 8; 


LM. (o), 


10. M, (z) < z 


11. end for; 
12. end for;// 计 算 每 个 像素 点 的 局 部 幅 值 特征 
13. fori & 1 to K 


"EOS 
oer 


E (M, (z)- Miu) 
(mxn)/ K 


14. end for // AES Ael lg fUr GE RO 1 Sj aea 


15.fori*-1to mxn; 


| Im(F., (2) 
Re(F,, (z)) 


0,, (z)) €- tan 
16. if 0«6,(2)«90]||180«6,, (z) «270; 
17. Pio (z)€0; 


18. endif; 


19. if 90«6,, (z) «180 || 270 « &,, (z;)<360 ; 


20. P (u)cl; 

21. endif; 

22. forjc€lto 8; 

23. FG) € B, G)* (5G) XORE, (2;)2^ ; 
24. end for; 


25. end for; /计算 每 个 像素 点 的 局 部 相位 特征 
26. fori & 1 to K; 


s e Dda (©) 
(mxn)/K 


D E 


(mxn)/ K 


xd = 
Ln 


27. end for; // 计 算 每 个 子 块 的 相位 特征 的 均值 与 标准 差 
28. fori & 1 to K; 


29. MIX, ,; -W, (M. Mr ) +W, (Pfa > D ) ; 


30. end for; / 幅 值 特征 与 相位 特征 的 加 权 融 合 
4 ”仿真 及 实验 结果 分 析 


4.1 实验 环境 和 参数 设置 

实验 需要 选取 合理 的 图 像 分 块 数 大 小 以 及 本 文 算法 融合 
特征 的 权 值 参数 大 小 。 实 验 的 评价 指标 为 平均 人 脸 识 别 率 以 
及 特征 提取 平均 时 间 。 将 本 文 算法 与 基于 Gabor 幅 值 与 相位 
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的 变化 。 


三 种 人 脸 数据 用 


图 像 示 例如 图 5 所 示 。 


Tr 


三 种 人 脸 数 据 集 图 像 示例 


Fig.5 Examples of three face data set images 
实验 选取 Gabor 滤波 器 参数 为 c=2r ，kw T, f-42, 


2 

并 提取 图 像 的 五 个 尺度 和 八 个 方向 上 的 特征 信息 。 实 验 的 软 
硬件 环境 为 ， CPU 为 Intel? Core™ i3-2310M CPUG2.10 
Hz; 内 存 为 6.00 GB; OS: Windows 7 Ultimate 64 bit; SE 
验 的 仿真 工具 : MATLAB 2014(a). 
4.2 分 块 数 大 小 的 选取 
针对 基于 Gabor 小 波 变换 的 特 行 
大 小 对 于 算法 提取 的 特征 的 表征 能 力 影 响 较 大 。 因 为 分 块 数 
较 小 时 ， 子 块 图 像 较 大 ， 则 提取 出 的 特征 不 够 丰富 ， 而 当 分 
块 数 较 大 时 ， 子 图 像 较 小 ， 则 特征 维 数 较 高 ， 并 且 容 易 受 噪 
声 影 响 。 因 此 ， 针 对 三 种 不 同 的 人 脸 数据 集 ， 分 别 通过 本 文 
算法 、 基 于 Gabor 幅 值 与 相位 的 特征 提取 算法 、 基 于 Gabor 
幅 值 的 特征 提取 算法 进行 特征 提取 ， 并 采用 相同 的 主 成 分 分 
析 法 对 提取 的 特征 进行 降 维 ， mdb i 
行 识别 分 类 ， 计 算 不 同 分 块 大 小 情况 1 丛 识 别 
留 一 交叉 验证 法 , 对 不 同 分 坎 大 小 情况 下 的 每 种 算法 进行 
次 测试 , 计算 20 次 测试 的 平均 人 脸 识 别 率 。 图 6 为 不 同 分 块 
大 小 情况 下 的 人 脸 识别 率 的 比较 。 

从 图 6 中 可 以 看 出 ， 当 达到 某 一 分 块 数 后 ， 本 文 算法 与 
基于 Gabor 幅 值 与 相位 的 特征 提取 算法 使 得 后 续 的 人 脸 识别 
率 持续 达到 最 高 且 没 有 下 降 趋势 ， 而 基于 Gabor 幅 值 的 特征 
提取 算法 的 后 续 人 脸 识 别 率 波 动 较 大 。 这 主要 是 因为 对 于 光 
照度 较 差 的 人 脸 图 像 而 言 ， 诸 如 眼睛 、 鼻 子 、 嘴 巴 等 脸 部 特 
征 的 影响 并 不 明显 。 当 分 块 数 较 大 时 ， 虽 然 其 连贯 性 受到 一 
定 程度 的 破坏 ， 但 Gabor 相位 信息 对 光照 的 补偿 能 力 很 强 ， 
所 以 本 文 算 法 与 基于 Gabor 幅 值 与 相位 的 特征 提取 算法 的 识 


图 5 


FE 提 取 算 法 ， 合 适 的 分 块 


ChinaXiv 合 作 期 刊 


录用 定稿 AER, Ai 改进 的 Gabor 小 波 变换 特征 提取 算法 第 37 卷 第 2 期 

的 特征 提取 算法 023、 基 于 Gabor 幅 值 的 特征 提取 算法 外 作对 征 提 取 算 法 、 基 于 Gabor 幅 值 与 相位 的 特征 提取 算法 和 本 文 

比 ， 比 较 在 光照 和 姿态 等 外 界 因 素 干 扰 情况 下 算法 提取 的 特 ”算法 的 最 佳 分 块 数 分 别 为 365、32、16; 在 Yale B 数据 集 下 ， 

征 对 于 后 续 人 脸 识别 率 的 影响 ， 其 次 比较 三 种 算法 的 特征 提 三 种 算法 的 最 佳 分 块 数 分 别 为 64. 49. 32; 在 CMU-PIE 数 

取 的 平均 时 间 。 实 验 选取 Yale. Yale B. CMU-PIE 三 种 人 脸 。” 据 集 下 ， 三 种 算法 的 最 佳 分 块 数 分 别 为 64、49、25。 

数据 集 来 比较 分 析 三 种 算法 。Yale 人 脸 数 据 集 包 含 了 15 个 100 

人 165 张 人 脸 图 像 ， 其 中 包括 了 光照 、 表 情 和 姿态 的 变化 ; 98 Lacon 

Yale B 人 脸 数 据 集 包 含 了 不 同人 种 的 38 个 人 2 432 张 人 脸 图 ~ o% 本 
H ES Gabor 幅 值 + 相位 

像 ， 其 中 包括 光照 与 姿态 的 变化 ; CMU-PIE 人 脸 数据 集 包 含 © 94 ——— A 

T 68 个 人 3332 张 人 脸 图 像 ， 其 中 包括 了 光照 、 表 情 和 姿态 x 
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(a)Yale 数据 集 
100 
98 
^. 96 
x% 
2 —Á— Gabor 幅 值 
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>a 一 一 本 文 算法 
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DI VS MOS 
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100 
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— 96 
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né 92 一 一 Gabor 幅 值 + 相位 
50 一 * 一 本 文 算法 
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分 块 个 数 CIO 
(c)CMU-PIE 数据 集 
图 6， 不 同 分 块 大 小 情况 下 的 人 脸 识 别 率 的 比较 


Fig.6 Comparison of face recognition rate under different block size 


4.3 ”特征 融合 的 权重 参数 选取 


验 设置 不 同 
3.3 节 可 知 ， 
者 和 为 1， 


数据 集 ， 分 别 通过 设置 了 不 同 特征 


别 率 没 有 受到 影响 。 从 图 中 还 可 以 看 出 ， 本 文 算法 与 基于 
Gabor 幅 值 与 相位 的 特征 提取 算法 相 比 ， 较 小 的 分 块 数 就 能 
使 人 脸 识别 率 达 到 最 高 ， 这 主要 是 因为 本 文 算法 考虑 了 幅 值 
特征 的 局 部 关联 性 。 从 图 6(a) 中 可 以 看 出 ， 使 得 识别 率 达 到 
最 高 时 , 基于 Gabor 幅 值 的 特征 提取 算法 的 最 佳 分 块 数 为 36 
个 ， 基 于 Gabor 幅 值 与 相位 的 特征 提取 算法 的 分 块 数 为 
32-128 个 ， 而 本 文 算法 的 分 块 数 为 16~128 个 。 因 为 分 块 数 
越 多 ， 时 间 代 价 就 越 大 ， 所 以 三 种 算法 的 最 佳 分 块 数 均 取 最 
小 值 。 同 理 可 以 得 到 图 6(b) 中 本 文 算法 、 基 于 Gabor 幅 值 与 
相位 的 特征 提取 算法 、 基 于 幅 值 的 特征 提取 算法 的 最 佳 分 块 
数 分 别 为 32、49、64，(c) 中 三 种 算法 的 最 佳 分 块 数 分 别 为 
25、49、64。 综 上 ， 在 Yale 数据 集 下 ， 基 于 Gabor 幅 值 的 特 


Ez 


的 特征 权重 


E, JB 
幅 值 特征 权重 为 Wy , 
特征 权重 参数 的 范围 为 0. 


由 于 本 文 算法 需要 融合 Gabor 幅 值 与 相位 特征 ， 
lioe E doc SUR 


1-09, 


征 提 取 j 


X) 


JH 


输入 到 相 
别 率 。 采 用 


不 同 特 生 


ZE ABAH XT Ai 
化 的 因素 ， 


留 一 


低 。 


Ly 当 幅 


同 的 支持 向 量 机 中 进行 识 另 
交叉 验证 法 ， 
进行 20 次 测试 ， 计 算 20 次 测试 的 平 ] 
FE 权重 参数 下 的 人 脸 识别 率 比 

从 表 1 可 以 看 出 
相位 特征 的 权重 参数 W, 偏 小 


值 特 征 


对 每 一 


所 以 实 


参数 。 根 据 
相位 特征 权重 为 W,， 且 两 
将 三 


种 不 同 的 人 


权重 的 本 文 算法 进行 特 


I 分 类 ， 
种 特征 


较 。 


同 的 主 成 分 分 析 法 对 提取 的 特征 进行 降 维 ， 
计算 得 到 人 脸 识 
权重 参数 的 设置 
名人 脸 识 别 率 。 表 1 为 


的 权重 参 


数 W, 1 
时 , 对 于 三 种 数据 集 的 人 脸 识 别 


局 大 ， 而 


主要 是 区 


以 需要 相位 4 


RE IE 


而 相位 特征 


DE: 


对 


为 这 三 种 数据 集 
于 光照 变化 


EC, 偏 大 


tm M4 
Lis —À 


在 Yale 数据 集 . 


W,-03H 


W, =0.7 RF Wy -04 H. 
上 的 识别 率 达到 最 高 


i, AR 


中 均 包 含 了 光照 变 
有 较 好 的 鲁 棒 性 ， 所 


中 还 可 以 看 


wW, =0.6 时 ， 本文 算法 
; 4 W, =0.4 H W, 2061, 


AEG, F: 


本 文 算法 在 Yale B 数据 集 上 的 识别 率 达 到 最 高 ; 24 Wy -03 H 


W, =07 或 者 W, =0.4 E] W, 
文 算法 在 CMU-PIE 数据 集 上 的 识别 率 达 到 最 高 。 为 了 统 
起 见 ， 选 定 幅 值 特征 的 权重 参数 W, 为 0.3， 相 位 特征 的 权 寻 
参数 W 为 0.7。 


0.6 或 者 W, =05 且 W=05 时 ， 本 


| 


my 


改进 的 Gabor 小 波 变换 特征 提取 算法 


ChinaXiv 合 作 期 刊 
第 37 卷 第 2 期 


BYA 


加 
以 


ra 
R 


均 时 间 是 高 于 基于 Gabor 幅 值 的 特征 提取 算法 的 ， 
了 相位 特征 的 提取 时 间 。 从 以 上 所 有 的 实验 结果 可 以 得 出 
下 结论 ， 本 文 算法 在 总 体 表现 上 是 优 于 其 他 两 种 算法 的 。 
然 本 文 算法 的 特征 提取 平均 时 间 高 于 基于 Gabor 幅 值 与 相 


表 1 不 同 特征 权重 参数 下 的 人 脸 识别 率 比 较 
Table ] Comparison of face recognition rate under different feature 
weight parameters 
特征 权重 人 脸 识别 率 /% 
W, W, Yale 数据 集 Yale B 数据 集 CMU-PIE 数据 集 
0.1 0.9 98.16 97.58 98.84 
0.2 0.8 98.62 98.12 99.07 
0.3 0.7 99.50 98.89 99.72 
0.4 0.6 99.50 99.56 99.72 
0.5 0.5 97.15 98.02 99.72 
0.6 0.4 96.72 96.97 98.16 
0.7 0.3 96.35 95.98 97.57 
0.8 0.2 95.52 95.21 97.34 
0.9 0.1 95.17 95.04 96.20 


4. 算法 提取 特征 的 表征 能 力 的 比较 


实验 将 三 种 不 同 的 人 脸 数据 集 ， 分 别 通过 本 文 算法 、 基 


EN 


TH 


于 Gabor 由 


值 与 相位 的 特征 提取 算法 、 基 于 Gabor 幅 值 的 特 


下 提取 算法 进行 特征 提取 ， 并 采用 相同 的 主 成 分 分 析 法 对 提 
区 的 特征 进行 降 维 ， 输 入 到 相同 的 支持 向 量 机 中 进行 识别 分 
类 ， 计 算得 到 人 脸 识别 率 。 
法 进行 20 次 测试 ， 计 算 20 次 测试 的 平均 人 脸 识别 率 。 表 2 


采用 留 一 交叉 验证 法 ， 对 每 种 算 


为 不 同 算法 提取 特征 的 表征 能 力 的 比较 。 
表 2 不 同 算法 提取 特征 的 表征 能 力 的 比较 /% 
Table 2 Comparison of characterization capabilities for different 
algorithms/% 
数据 集 Gabor 幅 值 Gabor 幅 值 + 相位 本 文 算法 
Yale 95.06 98.25 99.50 
Yale B 96.31 98.62 99.56 
CMU-PIE 96.93 99.08 99.72 
从 表 2 可 以 看 出 ， 本 文 算法 提取 的 特征 具有 最 好 的 表征 


能 力 ， 在 三 种 算法 中 具有 最 高 的 人 脸 识别 率 。 本 文 算法 的 特 


M 


fi 


E 向 量 表征 能 力 优 于 基于 Gabor 
法 ， 主 要 是 因为 本 文 算法 考虑 了 幅 值 特 4 


i [E 与 相位 的 特征 提取 算 
的 局 部 关联 性 ， 而 


T 


本 文 算 法 的 加 权 融 合 特征 优 于 特征 的 级 联 融 合 。 基 于 


Gabor 幅 值 与 相位 的 特征 提取 算法 的 特征 向 量 表征 能 力 优 于 


基于 Gabor 


融 值 的 特征 提取 算法 ， 主 要 是 因为 相位 特征 对 于 


光照 的 变化 具有 很 好 的 鲁 棒 性 ， 所 以 相位 特征 的 融入 使 得 算 


法 的 鲁 棒 性 得 到 了 提升 。 综 上 所 述 ， 本 文 算法 提取 的 特征 向 


幅 值 与 相位 的 特征 提取 算法 、 
法 的 特征 提取 平均 时 间 。 计 算 每 一 种 算法 分 别 在 三 种 数据 集 
下 的 特征 提取 平均 时 间 。 三 种 算法 的 特征 提取 平均 时 间 比 较 
如 表 3 Bras. 


基于 Gabor 
法 在 基于 Gabor | 
值 特征 的 计算 以 及 加 权 融 合 特征 的 计算 ， 所 以 增加 了 计算 时 
间 。 显 然 基 于 Gabor 幅 值 与 相位 的 特征 提取 算法 的 特征 提取 


量 的 表征 能 力 是 优 于 其 他 两 种 算法 的 。 
4.5 特征 提取 平均 时 间 的 比较 


本 节 实 验 在 三 种 数据 集 下 ， 比 较 本 文 算 法 、 基 于 Gabor 


基于 Gabor 幅 值 的 特征 提取 算 


a 


t 


从 表 3 可 以 看 出 ， 本 文 算法 的 特征 提取 平均 时 间 略 高 于 
局 值 与 相位 的 特征 提取 算法 。 主 要 是 因为 本 文 算 
值 与 相位 的 特征 提取 算法 中 加 入 了 局 部 幅 


位 的 特征 提取 算法 和 基 
在 最 终 的 人 脸 识 别 率 上 本 文 算法 高 于 这 两 种 算法 。 


在 


于 Gabor 幅 值 的 特征 提取 算法 ， 但 是 


保证 特征 提取 平均 时 间 略 有 增加 的 同时 ， 提 高 了 在 光照 、 


姿 


态 等 外 界 因素 干扰 情况 下 最 终 的 人 脸 识 别 率 ， 很 好 地 增强 


了 


特征 的 表征 能 力 ， 降 低 了 特征 向 量 的 维度 。 


表 3 三 种 算法 的 特征 提取 平均 时 间 比 较 
Table 3 Comparison of average time of feature extraction of three 
algorithms 
"T E 寺 征 提 bs l 
Gabor 幅 值 Gabor 幅 值 + 相位 。 本 文 算法 

Yale 数据 集 0.15 0.36 0.41 
Yale B 数据 集 0.11 0.29 0.38 
CMU-PIE 数据 集 0.10 0.25 0.35 


5 


结束 语 


为 了 解决 基于 Garbor 小 波幅 值 与 相位 的 人 脸 特 征 提取 


方法 的 级 联 特征 融合 方式 使 得 特征 向 量 维 度 较 高 的 问题 ， 本 


文 提出 了 一 种 改进 
算 局 部 幅 值 特征 与 局 部 相位 特征 来 增强 特征 之 间 的 关联 性 ; 
L5 


将 


本 


果 
特 
下 
的 


的 Gabor 小 波 变换 特征 提取 算法 。 通 过 计 


局 部 幅 值 特征 与 局 部 相位 特征 加 权 融 合 相 比 于 两 者 级 联 ， 
氏 了 提取 特征 的 维度 并 增强 了 特征 的 表征 能 力 。 从 实验 结 
来 看 ， 本 文 提 出 的 算法 所 提取 的 特征 有 较 高 的 表征 能 
征 维度 较 低 ， 但 是 特征 提取 的 平均 时 间 相对 较 高 。 因 此 ， 
一 步 工作 将 会 考虑 降低 算法 的 时 间 复 杂 度 ， 减 少 特征 提取 
平均 时 间 。 
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